Search for: All records

Abstract The marine atmospheric boundary layer (ABL) and oceanic boundary layer (OBL) are a two-way coupled system. At the ocean surface, the ABL and OBL share surface fluxes of momentum and buoyancy that incorporate variations in sea surface temperature (SST) and currents. To investigate the interactions, a coupled ABL–OBL large-eddy simulation (LES) code is developed and exercised over a range of atmospheric stability. At each time step, the coupling algorithm passes oceanic currents and SST to the atmospheric LES, which in turn computes surface momentum, temperature, and humidity fluxes driving the oceanic LES. Equations for each medium are time advanced using the same time step but utilize different grid resolutions: the horizontal grid resolution in the ocean is approximately four times finer, e.g., (Δx_o, Δx_a) = (1.22, 4.88) m. Interpolation and anterpolation (its adjoint) routines connect the atmosphere and ocean surface layers. In the simplest setup of a statistically horizontally homogeneous flow, the largest scale ABL turbulent shear-convective rolls leave an imprint on the OBL currents in the upper layers. This result is shown by comparing simulations that use coupling rules that are applied either instantaneously at everyx–ygrid point or averaged across anx–yplane. The spanwise scale of the ABL turbulence is ∼1000 m, while the depth of the OBL is ∼20 m. In these homogeneous, fully coupled cases, the large-scale spatially intermittent turbulent structures in the ABL modulate SST, currents, and the connecting momentum and buoyancy fluxes, but the mean profiles in each medium are only slightly different. 

Recent high-resolution large-eddy simulations (LES) of a stable atmospheric boundary layer (SBL) with mesh sizes N=(5123,10243,20483) or mesh spacings ▵=(0.78,0.39,0.2) m are analyzed. The LES solutions are judged to be converged based on the good collapse of vertical profiles of mean winds, temperature, and low-order turbulence moments, i.e., fluxes and variances, with increasing N. The largest discrepancy is in the stably stratified region above the low-level jet. Subfilter-scale (SFS) motions are extracted from the LES with N=20483 and are compared to sonic anemometer fields from the horizontal array turbulence study (HATS) and its sequel over the ocean (OHATS). The results from the simulation and observations are compared using the dimensionless resolution ratio Λw/▵f where ▵f is the filter width and Λw is a characteristic scale of the energy-containing eddies in vertical velocity. The SFS motions from the observations and LES span the ranges 0.1<Λw/▵f<20 and are in good agreement. The small, medium, and large range of Λw/▵f correspond to Reynolds-averaged Navier–Stokes (RANS), the gray zone (a.k.a. “Terra Incognita”), and fine-resolution LES. The gray zone cuts across the peak in the energy spectrum and then flux parameterizations need to be adaptive and account for partially resolved flux but also “stochastic” flux fluctuations that represent the turbulent correlation between the fluctuating rate of strain and SFS flux tensors. LES data with mesh 20483 will be made available to the research community through the web and tools provided by the Johns Hopkins University Turbulence Database. 

Abstract A persistent spatial organization of eddies is identified in the lowest portion of the stably stratified planetary boundary layer. The analysis uses flow realizations from published large-eddy simulations (Sullivan et al. in J Atmos Sci 73(4):1815–1840, 2016) ranging in stability from near-neutral to almost z-less stratification. The coherent turbulent structure is well approximated as a series of uniform momentum zones (UMZs) and uniform temperature zones (UTZs) separated by thin layers of intense gradients that are significantly greater than the mean. This pattern yields stairstep-like instantaneous flow profiles whose shape is distinct from the mean profiles that emerge from long-term averaging. However, the scaling of the stairstep organization is closely related to the resulting mean profiles. The differences in velocity and temperature across the thin gradient layers remain proportional to the surface momentum and heat flux conditions regardless of stratification. The vertical thickness of UMZs and UTZs is proportional to height above the surface for near-neutral and weak stratification, but becomes thinner and less dependent on height as the stability increases. Deviations from the logarithmic mean profiles for velocity and temperature observed under neutral conditions are therefore predominately due to the reduction in eddy size with increasing stratification, which is empirically captured by existing Monin–Obukhov similarity relations for momentum and heat. The zone properties are additionally used to explain trends in the turbulent Prandtl number, thus providing a connection between the eddy organization, mean profiles, and turbulent diffusivity in stably stratified conditions. 

Ocean fronts are an important submesoscale feature, yet frontogenesis theory often neglects turbulence – even parameterized turbulence – leaving theory lacking in comparison with observations and models. A perturbation analysis is used to include the effects of eddy viscosity and diffusivity as a first-order correction to existing strain-induced inviscid, adiabatic frontogenesis theory. A modified solution is obtained by using potential vorticity and surface conditions to quantify turbulent fluxes. It is found that horizontal viscosity and diffusivity tend to be readily frontolytic – reducing frontal tendency to negative values under weakly non-conservative perturbations and opposing or reversing front sharpening, whereas vertical viscosity and diffusivity tend to only weaken frontogenesis by slowing the rate of sharpening of the front even under strong perturbations. During late frontogenesis, vertical diffusivity enhances the rate of frontogenesis, although perturbation theory may be inaccurate at this stage. Surface quasi-geostrophic theory – neglecting all injected interior potential vorticity – is able to describe the first-order correction to the along-front velocity and ageostrophic overturning circulation in most cases. Furthermore, local conditions near the front maximum are sufficient to reconstruct the modified solution of both these fields. 

A common technique for estimating the sea surface generation functions of spray and aerosols is the so-called flux–profile method, where fixed-height concentration measurements are used to infer fluxes at the surface by assuming a form of the concentration profile. At its simplest, this method assumes a balance between spray emission and deposition, and under these conditions the concentration profile follows a power-law shape. It is the purpose of this work to evaluate the influence of waves on this power-law theory, as well as investigate its applicability over a range of droplet sizes. Large-eddy simulations combined with Lagrangian droplet tracking are used to resolve the turbulent transport of spray droplets over moving, monochromatic waves at the lower surface. The wave age and the droplet diameter are varied, and it is found that droplets are highly influenced both by their inertia (i.e., their inability to travel exactly with fluid streamlines) and the wave-induced turbulence. Deviations of the vertical concentration profiles from the power-law theory are found at all wave ages and for large droplets. The dynamics of droplets within the wave boundary layer alter their net vertical fluxes, and as a result, estimates of surface emission based on the flux–profile method can yield significant errors. In practice, the resulting implication is that the flux–profile method may unsuitable for large droplets, and the combined effect of inertia and wave-induced turbulence is responsible for the continued spread in their surface source estimates. 

Abstract. The science guiding the EUREC4A campaign and its measurements is presented. EUREC4A comprised roughly 5 weeks of measurements in the downstream winter trades of the North Atlantic – eastward and southeastward of Barbados. Through its ability to characterize processes operating across a wide range of scales, EUREC4A marked a turning point in our ability to observationally study factors influencing clouds in the trades, how they will respond to warming, and their link to other components of the earth system, such as upper-ocean processes or the life cycle of particulate matter. This characterization was made possible by thousands (2500) of sondes distributed to measure circulations on meso- (200 km) and larger (500 km) scales, roughly 400 h of flight time by four heavily instrumented research aircraft; four global-class research vessels; an advanced ground-based cloud observatory; scores of autonomous observing platforms operating in the upper ocean (nearly 10 000 profiles), lower atmosphere (continuous profiling), and along the air–sea interface; a network of water stable isotopologue measurements; targeted tasking of satellite remote sensing; and modeling with a new generation of weather and climate models. In addition to providing an outline of the novel measurements and their composition into a unified and coordinated campaign, the six distinct scientific facets that EUREC4A explored – from North Brazil Current rings to turbulence-induced clustering of cloud droplets and its influence on warm-rain formation – are presented along with an overview of EUREC4A's outreach activities, environmental impact, and guidelines for scientific practice. Track data for all platforms are standardized and accessible at https://doi.org/10.25326/165 (Stevens, 2021), and a film documenting the campaign is provided as a video supplement.